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We investigate an integrated optical chip immersed in atomic vapor providing several waveguide 
geometries for spectroscopy applications. The narrow-band transmission through a silicon nitride 
waveguide and interferometer is altered when the guided light is coupled to a vapor of rubidium 
atoms via the evanescent tail of the waveguide mode. We use grating couplers to couple between 
the waveguide mode and the radiating wave, which allow for addressing arbitrary coupling positions 
on the chip surface. The evanescent atom-light interaction can be numerically simulated and shows 
excellent agreement with our experimental data. This work demonstrates a next step towards 
miniaturization and integration of alkali atom spectroscopy and provides a platform for further 
fundamental studies of complex waveguide structures. 


Over the past decades, alkali atoms were not only sub¬ 
ject to a broad spectrum of research fields but also found 
their way into technological applications. Unlike solid 
state systems, the dispersion free properties of atoms 
are ideal for sensing and referencing tasks. While ultra 
cold atomic gases are best suited for ultra precise mea¬ 
surements and atomic clocks w, they usually require 
a large apparatus to do so. In contrast, devices based 
on thermal atomic vapors offer less precision, but have 
been successfully miniaturized and integrated for appli¬ 
cations e.g. in magnetometry EH3, frequency referencing 
mm, or atomic clocks m- Several approaches exist to 
achieve atom-light interaction on a microscopic scale such 
as atoms in hollow core fibers mm, nanofibers mm, 
micro- and nanocells mm, or anti-resonant reflecting 
optical waveguides (ARROW) on a chip |T8j [19]. Re¬ 
cently, an important step towards integration was made 
by surrounding solid core optical waveguides on a chip by 
a rubidium vapor cladding [20] , where the evanescent tail 
of the light mode interacts with the atomic vapor in close 
vicinity to the waveguide. With the existing technology 
of photonic integrated circuits on an industrial scale, this 
approach offers ideal conditions for using atomic vapors 
in e.g. sensing and communication applications. Light 
sources, interconnections, photonic devices and detectors 
can all be contained in a single chip, allowing for complex 
network and multiplexing designs, potentially even at the 
single photon level. Due to the small mode area of the 
evanescent field efficient atom light coupling is achieved 
and saturation can be reached already at low power. In 
this letter we increase the level of complexity compared 
to previous works m by adding Bragg couplers, curved 
waveguides and beam splitters to our silicon nitride pho¬ 
tonic structures. We show evanescent atom-light cou¬ 
pling by means of a simple waveguide strip as well as a 
Mach Zehnder type interferometer (MZI), providing both 
transmission and phase shift information. Our experi- 



FIG. 1. Experimental setup and layout of the structures, (a) 
Schematic of the experimental setup. For details see main 
text, (b) Microscopy image of the waveguide and (c) MZI 
structures, (d) Sketch of the MZI design and layer composi¬ 
tion. The coverage of the front coupler is depicted sliced to 
reveal the grating coupler. 


mental data can be well described by a theoretical model 
stemming from total internal reflection spectroscopy m 

The substrate of our optical chip consists of a 4 mm 
thick 1.5 inch diameter fused silica vacuum window, cov¬ 
ered with a 180 nm thick layer of silicon nitride (S^N^. 
The photonic structures are created in this layer by 
electron-beam lithography and subsequent dry etching. 
Details on the fabrication process can be found in our 
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previous paper [22]. Focusing grating couplers are used 
for in- and out-coupling of light. This type of coupler can 
be placed arbitrarily on the two dimensional chip surface, 
therefore allowing a larger number of individual devices 
on a single chip and more flexibility as compared to e.g. 
butt coupling from the side of the chip where only a one 
dimensional distribution is possible. Ah devices are com¬ 
pletely covered with a silicon oxide (SiCU) layer, except 
for the regions where we want the atoms to interact with 
the light held. Additionally, we deposit a 100 nm thick 
opaque layer of aluminum on top of the SiC >2 layer above 
each grating coupler to avoid leakage of light through the 
coupler and the detection of fluorescence light from the 
atoms inside the chamber volume. This layer also acts 
as a local mirror and therefore increases the coupling ef¬ 
ficiency. An overview of the layer composition is shown 
in figure [ljd). The chip is mounted into a custom made 
CF flange via a metal seal (Helicohex) and connected to 
a UHV chamber with the structures facing the inside of 
the chamber. After pumping and baking the chamber 
to a pressure < 10 -8 mbar, a rubidium ampule inside a 
bellow connected to the chamber is broken. We control 
the rubidium vapor pressure by the temperature of the 
bellow (reservoir), whereas we keep the chamber temper¬ 
ature always at a higher level (typically AT = 20 K) to 
avoid condensation on the chip surface. 

In our experiments, we focus a 780 nm laser through 
the substrate onto the input grating coupler of a specific 
device (see figure [IJa)). The output coupler of this device 
is imaged onto a 100 fi m pinhole to suppress any back¬ 
ground light which is coming from the vicinity of this 
port. After the pinhole we detect the signal with a photo 
multiplier tube (PMT). For the chip used here we uti¬ 
lized standard couplers which are not yet optimized for 
mode matching between a focused Gaussian beam and 
the waveguide mode but rather for direct coupling from 
a fiber tip. With direct coupling a transmission of up to 
6 % is achieved for the presented structures, whereas a 
record coupling efficiency of -0.62 dB has been reported 
for grating couplers [23]. Using a focused Gaussian beam 
the overall efficiency is 2.5 x 10“ 3 of the light which is 
sent to the chip and detected at the PMT for the simple 
waveguide. The waveguide quality was measured in the 
NIR and a propagation loss as low as 21 dB/m was ob¬ 
tained [22] . In the visible regime the losses are usually 
slightly increased, however, they are still negligibly small 
especially compared to the absorption due to the atomic 
vapor. 

For a first characterization of our system and to ex¬ 
amine solely its absorptive properties, we investigated a 
simple waveguide as shown in figure [ljb) with a width 
of 1.1 /im and a height of 180 nm. By numerical sim¬ 
ulation of the mode profile for this geometry, we infer 
that approximately 17% of the TE mode interact with 
the atomic vapor, which is the preferably guided mode 
for our waveguide design. The uncovered length of the 
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FIG. 2. Absorptive features for the simple waveguide, (a) 
Absorption spectra of the Rb D 2 line for various atomic den¬ 
sities n g = 1.8 x 10 12 cm -3 (blue), 6.8 x 10 12 cm -3 (green), 
1.3 x 10 13 cm -3 (red). The traces are normalized to the off- 
resonant transmission. The lines in darker color show the 
fit of the theory to the data, (b) Optical depth as a func¬ 
tion of reservoir temperature and detuning calculated from 
our model. Dashed lines indicate the positions of the data 
from (a). At a temperature of 122°C an optical depth of 1 is 
achieved for the Rb 85 5 Si/ 2 , F= 3 —)> 5P 3 / 2 transition. 


waveguide which is exposed to the atoms is ^1.2 mm. 
In figure [2| a) we show the normalized transmission spec¬ 
trum of such a device while scanning over the Rb D 2 line 
for different atomic densities. For a reservoir temper¬ 
ature of 113°C we already achieve an optical density of 
0.5. In principle, higher densities can be achieved by sim¬ 
ply increasing the reservoir temperature. However, as de¬ 
scribed later, we observed additional losses in the devices, 
possibly due to some deposition of Rb on the waveguide 
surface. Therefore we kept the reservoir temperature well 
below the chip temperature to reduce the amount of con¬ 
densation. The line shape of the spectrum exhibits some 
distinct deviation from the well known conventional Rb 
D 2 spectrum. This is caused by the enhanced Doppler 
broadening due to the 1.6 times larger wave vector in the 
waveguide compared to free space propagation, and by 
the limited transit time of the atoms traveling through 
the evanescent field. 
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FIG. 3. Transmission through a waveguide MZI with 2.2 mm path difference, (a) The gray solid lines show the calculated 
transmission of the MZI without any coupling to the atoms with a free spectral range of FSR ~ 71 GHz. The insets contain 
the experimental traces (blue) for various atomic densities n g — 2.4 x 10 13 cm -3 (z), 1.4 x 10 13 cm -3 (zz), 0.6 x 10 13 cm -3 (zzz). 
The solid red line is the theory fit while the red dashed line is the calculated transmission including the atomic contributions 
outside the scan range of the laser. Far off resonance this curve smoothly approaches the empty MZI fringe, (b) Zoom in on 
the insets from (a) with the experimental data (blue) and the theory fit (red). 
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To model our experimental data, we start with a finite- 
element analysis (COMSOL) to obtain the mode profile 
and propagation constant of the waveguide. With this 
information we can calculate the effective susceptibility 
Xeff of the atoms surrounding the waveguide, as described 
in previous works mmm- In this calculation we ne¬ 
glect any saturation effects, but include Doppler broad¬ 
ening and transit time broadening as well as self broaden¬ 
ing. Next, we add a cladding material with the complex 
refractive index riRb = \J 1 + Xeff to the waveguide in 
the COMSOL simulation. By running a frequency sweep 
with the detuning A around the center of the D 2 line, 
we obtain the complex propagation constant /3 r b- The 
transmission spectrum of the waveguide can then be cal¬ 
culated as 


2.2 mm. The shorter arm (1) is completely covered with 
Si 02 , whereas the longer arm (2) is only partially cov¬ 
ered, therefore offering a 2 mm long region for the guided 
mode to interact with the atoms. The modes from both 
arms are combined with a second Y-branch and guided 
to an output coupler. Besides the phase difference due to 
the unequal arm lengths, the light in the uncovered arm 
is picking up some additional phase caused by the real 
part of the susceptibility of the surrounding atoms. Also 
some of the light in the uncovered arm gets absorbed 
on resonance due to the imaginary part of the suscepti¬ 
bility. In total, this leads to a dispersive modulation of 
the MZI transmission as shown in figure [3] for different 
atomic densities. 

The transmission of the MZI can be calculated as 


T = / 0 x e- 2Im(/3Rb)z ', (1) 

where Iq is the intensity of the in-coupled light, and L 
is the length of the interaction region. As shown in fig¬ 
ure [2|a), this model fits well to our experimental data, 
where the frequency scaling, the center of the detuning 
and the Rb density are free fit parameters. In figure [2jb) 
we use our model to extrapolate the optical depth (O.D.) 
for a range of temperatures. From this follows that we 
can reach an O.D. of 1 at 122°C for the Rb 85 5 Si/ 2 , 
F= 3 —>> 5P 3 /2 transition and a waveguide with a length 
of 1.2 mm. 

The second type of devices we investigated is a sub- 
mm Mach Zehnder interferometer as shown in figure [ljc) 
and (d). Here the in-coupled light is split by means of a 
50/50 Y-branch into two arms with a path difference of 


Tmzi 


^ 2e «(^i(G-h)+/3Rbh) 


( 2 ) 


which describes the interference at the combining Y- 
branch. Here Ui, l\ and U 2 , h are the light amplitudes 
and lengths of arm 1 and arm 2, respectively. The length 
of the interaction region is denoted by lj. The propaga¬ 
tion constant for a waveguide with Si 02 cladding, as it 
is the case for arm 1, is denoted with /?i, whereas /3 r b 
is the complex propagation constant for a Rb cladding, 
as described earlier. With </> 0 we account for a phase off¬ 
set due to a temperature dependent change of the arm 
lengths. A fit of this model to our data is also plotted in 
figure [3^b) and shows excellent agreement with the Rb 
density, the amplitudes and the phase offset as the only 
free fit parameters. Figure [3ja) shows the MZI trans- 
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FIG. 4. Additional phase shift in the MZI. (a) The bright 
traces show the phase shift extracted from the data for atomic 
densities n g — 2.4 x 10 13 cm -3 (red), 1.4 x 10 13 cm -3 (green), 
0.6 X 10 13 cm -3 (blue). The dark curves are the correspond¬ 
ing calculated phase shifts for the parameters obtained from 
the fits in figure [3] (b) Calculated additional phase shift by 
the atoms modulo 2tt as a function of detuning and reservoir 
temperature. Dashed lines indicate the positions of the data 
from (a). 


mission for a larger detuning range with and without 
contribution from the atoms. From the fitted curves we 
deduce, that the amplitude in arm 2 is approximately ten 
times smaller than the amplitude in arm 1 and decreas¬ 
ing during the course of the experiment, thus causing a 
smaller visibility than expected from a 50/50 beam split¬ 
ter. We attribute this behavior to the larger length of 
arm 2 and therefore higher propagation losses on the one 
hand. On the other hand, it seems that some of the Rb 
atoms stick to the uncovered waveguide surface, increas¬ 
ing the losses in this arm additionally. A similar behav¬ 
ior was also found with the simple waveguide structures 
where we observed decreasing transmission over time. 

The most intriguing feature of a Mach Zehnder inter¬ 
ferometer is of course its ability to measure phase shifts. 
We now can extract this phase shift A ip due to the pres¬ 
ence of the atoms from our data transforming equation [2] 


and subtracting the phase shift of the bare MZI: 

a _ -if Tmzi \Ui\ 2 ~ |^ r 2| 2 exp(-2Im(^R b )/i) \ 

^ C ° S ( 2U1U2 exp (— Im(^R b )ii) J 

— [Pih — (Pi(l 2 — h) + Mi)} , 

where /3 q is the propagation constant in the waveguide 
without cladding (vacuum). Figure [4ja) shows the phase 
shifts corresponding to the data in figure [3j For the data 
with the highest atomic density of n g = 2.4 x 10 13 cm -3 , 
the light experiences an additional phase shift of up to 
0.15 x 7T. In figure [4jb) the calculated atomic phase shift 
for this particular device is shown in dependence of the 
reservoir temperature. Again we can extrapolate from 
our model, that an additional phase of it is reached at a 
temperature of 160° C, corresponding to an atomic den¬ 
sity of 1.7 x 10 14 cm -3 . Naturally a higher density is 
accompanied with strong absorption, but since the abso¬ 
lute value of the real part of the susceptibility is largest at 
the wings of the absorption lines, the off-resonant phase 
shift is still present, without much attenuation. The 
phase sensitivity of these devices can easily be increased 
by lengthening the uncovered arm or by decreasing the 
mode confinement. 

In conclusion, we have presented a hybrid system con¬ 
sisting of thermal alkali vapor and integrated photonic 
structures on a chip. Our optical chip houses several de¬ 
vices featuring grating couplers for flexible addressing of 
the individual devices. In future chip designs we are going 
to match the design of these couplers to our experimental 
conditions to improve the coupling efficiency. The trans¬ 
mission spectra of the simple strip waveguide revealed ab¬ 
sorption of the evanescent tail for various atomic densities 
and showed line broadening due to the Doppler effect and 
the short transit time of the atoms through the evanes¬ 
cent field. To reach more narrow lines, two photon spec¬ 
troscopy (e.g. electromagnetically induced transparency, 
EIT) can be utilized to cancel the Doppler shift, whereas 
the transit time broadening could be reduced by adding 
buffer gas. In addition to absorptive measurements we 
performed phase sensitive measurements using an inte¬ 
grated Mach Zehnder interferometer and could extract 
the atomic phase shift from our data for different Rb den¬ 
sities. By numerically simulating the light propagation 
in a waveguide surrounded by an atomic vapor cladding 
with a complex effective refractive index, we could re¬ 
produce the experimental data for both types of devices 
with excellent agreement. Over time we witnessed some 
degradation of the waveguide structures, possibly due 
to some build up of a Rb layer on their surface. This 
process appears to be partially reversible: after we cool 
down the Rb reservoir but keep the chip at ^ 200° C, the 
transmission increases again. For future chip generations 
we will investigate alternative materials and protective 
coatings to increase the lifetime of the devices. Addi¬ 
tionally we aim to include the optical chips in a vapor 
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cell using anodic bonding [25], which is a further step 
towards integration and miniaturization and also allows 
for a better temperature control, therefore reducing the 
risk of Rb condensation on the waveguides. This hy¬ 
brid system opens the door for future experiments with 
various waveguide geometries. The small mode area of 
the evanescent field on the order of A along the waveg¬ 
uide allows for systematic studies of interaction effects 
like self broadening in the one dimensional case, as it has 
been shown for two dimensions in a thin vapor cell m- 
These interactions could also be utilized to add nonlin¬ 
earity e.g. to a system of coupled ring resonators, which 
on their own create a synthetic gauge field for photons in 
the non-interacting regime [26 . 
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